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ABSTRACT 

We present ~ 2" resolution CO (3-2), HCO + (4-3) and 880/^m continuum images of the luminous 
infrared galaxy NGC 6240 obtained at the Submillimeter Array. We find that the spatially resolved 
CO (3-2), HCO + (4-3) and the 880/xm emission peaks between the two nuclear components that are 
both known to harbor AGNs. Our Large Velocity Gradient (LVG) analysis performed on each velocity 
channel suggests that the peak of the molecular gas emission traced in our observations is warm (T = 
20 - 100 K), dense (nn 2 = 10 5 0-5 4 cm -3 ) and moderately optically thin (r = 0.2 - 2) in the central 

1 kpc. We also find large column densities of ~ 10 23 cm~ 2 . Such extreme conditions are observed 
over ~ 300 km s _1 centered around the CO derived systemic velocity. The derived molecular gas 
mass from the CO (3-2) emission and a CO-to-H2 conversion factor commonly used for ULIRGs is 
(6.9±1.7) x 10 9 Mq, and this is consistent with the mass derived from previous CO (2-1) observations. 
The gas is highly turbulent in the central kpc (Aufwzi ~ 1175 km s -1 ). Furthermore, possible inflow 
or outflow activity is suggested from the CO (3-2) velocity distribution. We tentatively state that 
3.5 x 10 8 Mq of isolated CO (3-2) emission seen west of the northern disk may be associated with 
outflows from starburst superwinds, but the gas outflow scenario from one of the central AGN is not 
completely ruled out. Piecing all of the information together, the central region of NGC 6240 harbors 

2 AGNs, ~ 10 10 Mq of molecular gas mass, 5 x 10 7 M Q of dust mass, and has possible evidence of 
inflow and outflow activity. 

Subject headings: galaxies: formation, galaxies: interaction, galaxies: starburst, galaxies: kinematics 
and dynamics, galaxies: individual (NGC 6240) 



1. INTRODUCTION 

(Ultra) Luminous Infrared Galaxies [(U)LIRGs] are 
systems discovered in the local universe that emit 
extremely large luminosit y in the far infrared (see 
Sanders fc Mirabel 119961 : lLonsdale. Farrah fc Smith! 
20061 for review). Their disturbed optic a l/NIR 
morphologies (e.g . lArmus. Heckman fc Milevl 119871: 



gas physics and the evolution of a central black hole 
have tied the merger event with radial gas inflow, 
central starbur sts, black hole feeding, and a p o werful 
gas outflow (Springcl, Di Mattco & Hernq uistl 120051: 
Hopkins et al.l 120051: [Di Matteo. Springel fc Hernquistl 



Surace et all [199 
20021) 



iFarrah et al.l l200ll iBushouse et all 
highly compact molecular gas emis- 
suggest that colli- 



and the 

(jDownes fc Solomon! Il998h 



sion 

sions/mergers of gas rich galaxies and subsequent radial 
inflow of gas toward the central region is primarily 
responsible for triggering extens ive bursts of nuclear star 
formation (IScoyille et al.l l2000h. feeding th e obsc ured 
central AGN (llmanishi, Dudley fc Malonevl |2006| ). or 
most likely, both (|Genzel et al.lll998f) . Such inflow activ- 
ity was already predicted in early numerical simulation 
studies t hat involve collisions of two massive gas rich 
galax i es dBarnes fc Hernquistl [19 96: Mi hos fc Hernquistl 
119961: llono. Yun fc Mihosl 120041 ). More recent sim- 
ulation studies that include a more realistic set of 
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2005: iNaravanan et al.ll2006f) . These simulations have 
predicted a self-consistent evolutionary scenario from 
mergers to ULIRGs with a brief optical QSO phase 
([Hopkins et alJl2006f ). 

High resolution observations of molecular gas tracers 
in (U)LIRGs and have been limited to low J transition 
CO emission obtained at mm interferometers. Past sin- 
gle dish measurements of dense gas tracers such as HCN 
(|Gao fc Solomonir2004h have found that molecular gas in 
(U)LIRGs is highly compressed and dense. Thus, ob- 
taining sensitive and high spatial resolution images of 
high excitation molecular gas emission will provide ex- 
tremely valuable information regarding the star forma- 
tion process during a merger and will further constrain 
the physics implemented in the evolution simulations. In 
addition, recent surveys of molecular gas traced in high 
J transition CO in high-z galaxies (jGreve et al.l 120051 : 
iTacconi et al.l l2006) will allow us to conduct a direct com- 
parison of the gas properties with local counterparts. 

Here we present interferometric CO (3-2), HCO + (4- 
3), and 880/im continuum observations of the luminous 
infrared galaxy (LIRG) NGC 6240 obtained at the SMA 
(|Ho. Moran fc Loll2004h . Both molecular lines are known 
to trace warmer and denser molecular gas than the com- 
monly observed mm lines such as CO (1-0) or CO (2-1). 
Past SMA observations toward infrared bright galaxies 
suggest that the distribution of the CO (3-2) emission 



2 



is si gnificantly different from the lower t ransition lines of 
CO (jlono et al.ll200i iWang et al.ll20?a . 

1.1. Multi-Wavelength Properties of NGC 6240 

NGC 6240 is an infrared luminou s merger (Lfir = 
3.5 x 10 11 L Q ) (lYun fc Carillil 120021) at a distance of 
D L = 107 Mpc (1" = 493 pc) 7 . Its disturbed opti- 
cal morphology and the compact double nuclei suggest 
that NGC 6240 is a late stage merger. The presence 
of AGN and extended starbursts were already suggested 
from early X-ray obser vations using a v ariety of space- 
borne instruments (i.e. Ilkebe et"aTll2000L and references 
therein). More recently, Chandra has reso lved the central 
regio n into two AGNs separated by ~ 1" (Komoss a et al.l 
2003), providing significant evidence that NGC 6240 is 
a merger of two AGN host galaxies. Recent XMM- 
Newton spectroscopy has revealed a significant contribu- 
tion of starburst activity to the extend ed X-ray emission 
(|Boller et al.ll2003t iNetzer et al.ll2005ft . 

Optical and infrared observations using ground based 
and spaceborne telescopes have also been used to study 
this galaxy intensively. HST and Keck imaging have 
resolved each nucleus into several d istinct components 
(|Gerssen et alj 2004; M ax et al"1l2005f) with an apparently 
dusty central region. A large velocity gradient in ion- 
ized gas has been id entified using optical spectroscopy 
(|Gerssen et al.ll2004D . Near IR imagi ng has also shown 
the presence of two Jf- band nuclei (jScoville et al.ll2000L 
iBogdanovic et al"1l2003h . Signatu res of starbursts a ppear 
to dominate the MIR spectrum (|Lutz et alj [20031. Mid 
IR spectroscopy using IRS on the Spitzer Space Telescope 
suggest a minor (~ 3 - 20% ) contribution of th e AGN to 
the bolometric luminosity (jArmus et al J 12006). but mid 
infrared spectra obtained at Keck suggest that the dom- 
inant po wer-source could b e explained by starbursts or 
an AGN (jEgami et al.l l2006). Optical spectroscopy using 
Subaru has revealed significant H2 emission 0'.'2 north of 
the southern nucleus, suggesting that the medium be 
tween the two nuclei is dominate d by starburst super 
winds from the southern nucleus (jOhvama et al.l 2000) 
Finally, a star formation rate of 61 ± 30 M fi yr " 1 is ob 



tained from a fit to the FIR SED (|Yun fe Cariil 



Kcnnicutt 



2002h 



and using the Schmidt law derived in 
to covert to the Star Formation Rate (SFR). 

Interfcromctric millimeter studies have iden- 
tified signific ant amounts of mole c ular gas in 
NGC 6240 dBrvant fc Scovilld fl99flt iTacconi et aD 
1999; Nakanis hTet al.ll2005h . The bulk of cold and dense 
gas peaks between the two optical/X-ray nuclei, close to 
the r egion where significant H2 emission ([Oh vama et al.l 
120001) is found (Figure 1). Millimeter continuum 
emission also peaks in the same region. Molecular gas 
kinematics traced in CO (2-1) emission show a large 
velocity gradient between the two nuclei, suggesting 
that ~ 10 9 M Q of molecular gas is in the process of 
accumulating in the new dynamical center between the 
two nuclei. 

MERLIN observations (|Beswick et al.ll200lD have de- 
tected both AGNs with highly disturbed H I gas in ab- 
sorption. Further, high res olution VLBA observations 
(jGallimore fc B cswick 2004) have detected two compact 
radio sources (Nl and SI in Figure [1]) and an additional 

7 Adopting H = 70 km s" 1 Mpc^ 1 , U M = 0.3, n A = 0.7 



component (S2) 072 north o f the southern nucleus (SI). 
iGallimore fc Beswickl (|2004ft attribute Nl and SI to com- 
pact radio emission associated with the AGNs, and S2 as 
a radio supernova or ma terial ejected from SI. An H2O 
mase r detected near SI (|Hagiwara. Diamond fc Mivoshil 
2003) is attributable to either a jet associated with 
the southern AGN, or with shocks caused by the im- 
mense collisio n . Fro m their low resolution VLA maps, 
iColbert et al.l (|1994T ) have found significant centimeter 
emission west of Nl. This emission is resolved into 
four distinct components, and they suggest shell material 
from galactic super-wind as an origin of the emission. 

This article is organized as follows. In §2, we present 
our observational method and the detailed calibration 
scheme adopted for the SMA observations. The CO (3- 
2) and HCO + (4-3) moment maps are presented in §3, 
along with the 880^m continuum image which was made 
by averaging the emission in the line free channels. §4 
focuses on the physical properties of gas. We provide a 
discussion of our results in §5, and conclude this paper 
in §6. 

2. OBSERVATIONS AND DATA REDUCTION 

NGC 6240 was observed on October 9 and 14, 2005 us- 
ing both the compact (unprojected baselines =16-69 
meters) and extended configurations (unprojected base- 
lines = 50 - 182 meters) of the SMA. These data were 
taken in excellent atmospheric conditions (T230 = 0.05 - 
0.06). The digital correlator was configured with 2048 
channels for the entire 2 GHz bandwidth for each side- 
band. The receivers were tuned to the redshifted CO (3- 
2) (z-Vest = 345.796 GHz) emission line in the lower side- 
band (LSB), which allows us to simultaneously obtain 
the HCO+(4-3) {v rest = 356.734 GHz) line in the up- 
per sideband (USB). The adopted tuning frequency was 
337.533 GHz (LSB) using v sys = 7339 km s" 1 . We used 
q(J2000) = 16 h 52 m 58.9 s and j )(J2000) = 2°24'04.00" 
(jLaurent-Muehleisen et al.ll 1997; as our phase reference 
center. Initial data calibration were carried out in the 
IDL based SMA calibration package MIR, where 3C279, 
3C454.3 and Uranus were used for bandpass calibration. 
Absolute flux calibration was performed using Uranus, 
and time dependent gain calibration was performed us- 
ing the 1.6 Jy QSO 1743-038 (14 degrees away from 
NGC 6240) and the 1.4 Jy QSO 1549+026 (16 degrees 
away from NGC 6240) as references. 

Synthesis imaging was performed in MIRIAD. Chan- 
nel maps with 25 km s _1 velocity resolution were made 
for the CO (3-2) line. We further varied the visibility 
weighting from natural (i.e. maximum sensitivity) which 
gave a synthesized beam size of 2'.'0 x l'/8 (position angle 
= 49 degrees; Figure [2] (left)) to uniform (i.e. maxi- 
mum angular resolution) which gave a synthesized beam 
size of I'll x Of! 8 (position angle = 49 degrees; Figure [2] 
(right)). Unity weights are assigned for each visibility 
point in natural weighting, giving the best S/N and the 
coarsest angular resolution. On the other hand, uni- 
form weighting gives the highest angular resolution at 
the cost of lower S/N. This is achieved by weighting each 
visibility point by the reciprocal of the number of data 
points found in a grid centered around each visibility 
point. The rms noise level was 22 mJy and 78 mJy for 
the natural and uniformly weighted maps respectively. 
In addition, in order to derive the peak coordinates of 
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the CO (3-2) emission, we have constructed an image 
using only the longest baseline data which gave a syn- 
thesized beam size of 0'.'9 x 0'.'5 (position angle = 41 de- 
grees) . For the natural weighted beam, the first sidelobes 
are significant at the 50% level ~ 11" north and south 
of the main synthesized beam. Sidelobes at 10% level 
are seen in various positions, the effect of which should 
be negligible with proper deconvolution. For the uni- 
formly weighted beam, the first significant (~ 30%) side- 
lobes are located ~ 5" north-south of the main beam. 
CLEANing was performed down to 1.5 times the rms 
noise level of each channel map. For the deconvolution 
algorithm, we chose CLEAN instead of the Maximum 
Entropy Method (MEM) since the emission is strong and 
relatively compact with both the natural and uniformly 
weighted beams. 

The HCO + (4-3) channel maps were made using natu- 
ral weighting, which gave a 2'.'0 x l'/7 (position angle = 42 
degrees) beam. We note that the HCO + (4-3) emission 
occurs near the edge of the bandpass, and these observa- 
tions may be missing flux at the low velocity end. Con- 
tinuum subtraction was performed by removing a linear 
baseline in the visibilities in both sidebands. The 880 /im 
continuum image was made by adding all of the line free 
channels in the bandpass of both sidebands. The rms 
noise level for the continuum map was 6 mJy. From the 
S/N of the maps, we estimate the astrometric accuracy 
of the CO (3-2) maps presented here to be < O'.'l, and 
- O'.'l for the HCO+(4-3) and 880 ^m maps. 

3. RESULTS 

The channel maps of CO (3-2) emission are provided 
in Figure [3l Figure [4] (left) shows the integrated inten- 
sity map of the natural weighted CO (3-2) emission, and 
Figure 2] (right) shows the uniformly weighted image of 
the central 6". Figure [5] show the velocity distribution 
with the same angular resolution as in Figure [4] The 
velocity distribution of the western complex (WC here- 
after; see below) alone is shown in Figure [U The inte- 
grated intensity maps are clipped at 1.5 times the rms in 
each channel, whereas the velocity distribution maps are 
clipped at 3 times the rms. The position velocity maps 
obtained along the two AGNs are shown in Figure[7] Fig- 
ure [8] shows the HCO + (4-3) and 880/im emission images 
of NGC 6240. A summary of the derived parameters is 
shown in Table 1. 

3.1. CO (3-2) Emission 
3.1.1. Distribution 

The size of the CO (3-2) emission region in the inte- 
grated intensity map of Figure|H(Ze/i) extends 4.4 kpc in 
both north-south and east-west directions. The coordi- 
nates of the peak determined from the highest resolution 
map obtained by just using data in the extended con- 
figuration are a(J2000) = 16 h 52 m 58.9 s and 5(J2000) = 
2°24'03'.'70. This is coincident with the radio compo- 
nent (S2) at a(J2000) = 16 h 52 m 58. 8994 s and <5(J2QOO) = 
2°24' 03"5925 (error = 0.004") (|Gallimore fc Beswickl 
2004) to within O'.'l which is comparable to the astro- 
metric accuracy of the SMA. While the inner region of 
the map shown in Figure |4] (left) is not spatially resolved 
with this map, there are several distinct ~ 1 kpc scale 
features that dominate the emission in the periphery; (1) 



the emission that extends to the east, (2) the emission 
that extends to the north, and (3) the isolated emission 
west of the main CO (3-2) concentration (WC). 

The ~ 2 kpc eastern extension is significant in the ve- 
locity range of 6964 to 7489 km s _1 (see Figured]), and 
appears to originate from the southeastern edge of the 
central emission region. On the other hand, the ~ 2 kpc 
northern extension is significant in the velocity range 
7264 to 7564 km s _1 , and appears to display a slight 
curve toward the east. The isolated component (WC) 
was detected ~ 9" (4.4 kpc) west of the main CO (3-2) 
peak, and this emission is significant in the velocity range 
7114 to 7714 km s _1 . No significant optical counterpart 
is associated with WC. In addition, none of the past 
interferome tric CO observations have identified a simi- 
lar feature dBrvant fc Scovilldll999l : iTacconi et~aDll999L 
iNakanishi et al.l 2005[ ). but radio continuum observations 
have detected significan t extension of syn chrotron emis- 
sion in the same region (jColbert et al.lll99l . 

The total integrated intensity from the full resolution 
map including WC is 2834±20 Jy km s _1 , recovering 90% 
of the total flux measure d using a single d ish telescope 
(3180 ± 460 Jy km s" 1 ; iGreve et al.ll200l . Using the 
CO-H 2 conversion factor appropria te for ULIRGs (i.e. 
X = 0.8±0.2 Mq (K km s" 1 )" 1 pc 2 : lDownes fc Solomon! 

HS98D, and assuming h' co(3 _ 2) = L co{1 _ 0) , L co{3 _ 2) = 

(8.6 ± 0.1) x 10 9 K km s" 1 pc 2 translates to M H2 = 
(6.9 ± 1.7) x 10 9 Mq. The peak of the CO (3-2) emis- 
sion gives a peak H2 column density of Nh 2 = (2.0 ± 
0.3) x 10 23 cur 2 using X = 0.8 M (K km s" 1 )" 1 pc 2 . 
The peak brightness temperature observed with the ~ 2" 
beam is 8.6 K. The integrated intensity of WC alone 
is (143 ± 20) Jy km s _1 , which translates to Mh 2 = 
(3.5 ± 1.0) x 10 s M . 

The higher resolution image afforded by adopting uni- 
form weighting shown in Figure [4] {right) resolves out 
the extended structure, revealing the compact CO (3- 
2) emission that dominates the central 1.5 kpc. The 
integrated intensity derived from this map is 1839 ± 
20 Jy km s^ 1 and this is 65% of the total intensity de- 
rived from the full resolution image. The integrated in- 
tensity translates to Mh 2 = (4.5 ± 1.1) x 10 9 M Q assum- 
ing X = 0.8 ± 0.2 M Q (K km s" 1 )" 1 pc 2 . The CO (3-2) 
emission arises primarily in a compact structure centered 
between the two nuclei and slightly extended in the di- 
rection defined by the line connecting the two. 

3.1.2. Kinematics 

The total CO (3-2) FWZI linewidth of 1175 km s" 1 
is extremely large (see Figure [3]). The large linewidth, 
particularly in the medium surrounding the two AGNs, 
suggests extremely turbulent gas. While it is possible 
that the bulk of the gas seen in Figure [5] [left) is dom- 
inated by turbulent gas stirred during the massive colli- 
sion, the northeast-southwest velocity gradient centered 
around the northern nucleus (Nl) may be suggestive of 
a disk rotation. The higher resolution image in Fig- 
ure [5] {right) also shows this compact, rotation feature. 
Rapid changes in velocity of Av ~ 100 km s _1 are seen 
along the northern extended feature that appears to con- 
nect to this rotating disk, possibly suggesting collimated 
gas streaming toward the direction close to the line of 
sight. A steep velocity gradient like this is suggestive of 
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gas inflow or outflow. The velocity gradient in the vicin- 
ity of the southern nucleus (SI) is less obvious, possibly 
due to the dominance of local non-circular motion, or a 
face-on orientation of a rotating disk. Although our cur- 
rent data does not have high enough angular resolution to 
distinguish the exact origin, if the southern nucleus has 
a face-on rotating disk, then the large shift in velocity of 
~ 250 km s^ 1 that begins from the south of SI and along 
the eastern extension may also suggest inflow or out- 
flow activity. The velocity distribution (Figuie\El(right)) 
shows that the systemic velocities of the two nuclei are 
offset by ~ 100 km s _1 . The emission in WC is sig- 
nificant in the velocity range 7114 to 7564 km s" 1 (see 
Figure [3]), and does not show any significant evidence of 
a rotation in the moment map (Figure [6]). An overall ve- 
locity shift from the northern part to the southern part 
of WC is evident. 

The position- velocity diagrams (PVD) sliced along the 
two AGNs (Nl and SI) and centered exactly halfway 
between the two AGNs are shown in Figure for dif- 
ferent resolution maps. Figure [7] (left) clearly shows 
that the linewidth near the two nuclei is extremely large, 
with emission slightly lopsided toward the redshifted ve- 
locities (v > km s _1 ). The PVD obtained with the 
uniformly weighted map resolves out the extended emis- 
sion in the velocity extrema, showing more details in the 
central ±200 km s _1 (Figure [7] (middle)). Similar to 
Figure [7] (left), the peak of the emission is near the sys- 
temic velocity of 7339 km s — 1 . In contrast, the PVD 
obtained (Figure [7| (right)) with the highest angular 
resolution map completely resolves the central emission 
into two distinct peaks separated by ~ 70 km s _1 and 
~ 0'.'2 (100 pc). This double peaked morphology in the 
PVD suggests overlapping molecular cloud complexes in 
the medium between the two AGNs. The fact that these 
two peaks are nearly spatially coexistent, but separated 
in velocity further strengthens the argument that the gas 
kinematics in the central region is highly turbulent and 
non-circular. 

3.2. HCO + (4-3) Emission 

The distribution of the HCO + (4-3) emission (Fig- 
ure^ (left)) is highly compact, and it is concentrated in 
the central region with a northeast-southwest elongation 
of 4" x 2" (2x1 kpc), with an additional slight exten- 
sion toward the southeast. The peak of the HCO + (4- 
3) emission is a(J2000) = 16 h 52 m 58.9 s and <5(J2000) = 
2°24'03'/8, and this is consistent with the peak of the 
CO (3-2) to within the astrometric accuracy of the SMA. 
The integrated intensity is (28±4) Jy km s _1 , recovering 
40% of the total flux mea sured using a sing le dish tele- 
scope (78 ± 16 Jy km s' MGreve et all 120061 ) . The peak 
brightness temperature observed with the ~ 2" beam is 
0.5 K. The highly compact and low S/N emission pre- 
vents us from constructing an interpretable velocity dis- 
tribution map. The HCO + (1-0) emission obtained with 
a coarser synthesized beam also peaks at a similar lo- 
cation, although the emission is more extended in the 
northwest-southeast direction (Nakanishi et al. 2005). 

3.3. 880 fim Continuum Emission 

Similar to the CO (3-2) and HCO+(4-3) emission, the 
880/im emission also peaks between the two nuclei, near 



the radio continuum feature S2. The 880 fim peak is 
also consistent with the 1.3 mm peak within the uncer- 
tainties of both measurements (Figure [1]). The coordi- 
nates of the 880/xm peak are a(J2000) = 16 h 52 m 58.9 s 
and (5(J2000) = 2°24'03'/4. The extension to the east is 
coincident with the similar extension seen in the CO (3- 
2) map, but the continuum is much shorter. The size of 
the 880 nm emitting region is 6'/0 x 3'.'5 (2.9 x 1.7 kpc). 
The total flux detected from the SMA observations is 
(105 ±6) mJy, and this is consistent with the true single 
dish continuum measurement of 96 m Jy, which is calcu- 
lated from the SCUBA 850/rni flux (jKlaas et alJ l2001h 
and subtracting the contribution of the CO (3-2) emis - 
sion using the formulae given in ISeaauist et ail ()2004h . 
Assuming that all of the 880/im emission is associated 
with dust, the derived dust mass in the central 2 kpc 
is (4.7 ± 1.4) x 10 7 M© using the formula described in 
Wilson et al. (2006 in prep). This yields a gas-to-dust 
mass ratio of (147 ± 57) in the central 2 kpc region 
(see §3.1 for the gas mass). Although uncertainties in 
both mass calculations are large, this is in agreement 
with the values typically observed in the Galaxy. The 
peak of the continuum suggests that the bulk of the 
dusty starburst activity occurs in this region. This is 
in contrast to previous optical/NIR im aging studies (e.g. 
lPasquali~G allaghc r fc de G riis 2004) where the visible 
starburst activity is concentrated near the two optical/X- 
ray nuclei, suggesting that the SFRs (2.8 - 412 M Q yr" 1 ; 
iPasauali. Gallagher fc de G riis 2004) derived using these 
measurements are likely underestimated. 

4. ANALYSIS OF THE PHYSICAL PROPERTIES 

4.1. Large Velocity Gradient Modeling 

Our new CO (3-2) and HCO+(4-3) data allow us to 
investigate the physical conditions of the gas by taking 
the line ratios wit h previously pub l ished CO (1-0) and 
HCO+(1-0) data (jNakanishi et all l2005h . Past multi- 
line studies toward galaxies have been lim ited to line 
ratios obtained from single dish data (e.g. | Aalto et al.l 
119951: iPetitpas fc Wilsonll2000HYao et al.ll2003ft or the ra- 
tios between the integrated intensities for near by systems 
(e.g. iMatsushita et all 119991 : iMao et al.l [2000h . Here we 
present line ratios obtained in each velocity channel, and 
investigate the gas properties by comparing the ratios 
with the theoretical predictions obtained from Large Ve- 
locity Gradient (LVG) models dGoldreich fc Kwanill97i 
iTakakuwa. Mikami fc Saitolll998h . We" assume here that 
both CO and HCO + emission in the central kpc arise 
from the same molecular cloud. Although it is sug- 
gested that the filling factors of denser gas tracers may be 
smaller than the commonly observed CO (1-0) line (e.g. 
lAalto et al.lll995[ ). here we assume for simplicity that the 
(beam and volume) filling factors between CO (3-2) and 
CO (1-0), as well as HCO+(4-3) and HCO+(1-0) are 
similar. We note that, although the HCO + (4-3) emis- 
sion recovers only 40% of the total single dish flux (§3.2), 
the main focus of this analysis is to investigate the gas 
properties in the central ~ 1" where the missing flux 
should be negligible. By comparing the line ratios be- 
tween CO (3-2) and CO (1-0), and also HCO+(4-3) and 
HCO + (1-0), these assumptions allow a direct compari- 
son with LVG models that does not depend on the filling 
factor. 
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All data were convolved to the angular resolution of 
the CO (1-0) image (2'.'2 x 2'.'1; lowest angular resolution 
among the four data cubes), regridded to the same ve- 
locity resolution of 56 km s _1 , and two line ratio cubes, 
CO (3-2)/CO (1-0) (Rco) and HCO + (4-3)/HCO+(l- 
0) (Rhco+) were formed. The pixel-to-pixel correla- 
tion between the two ratio cubes were obtained and are 
shown as crosses in Figure [9] In order to show the re- 
gions with significant detection in all four transitions, 
only the data above two sigma in each channel are plot- 
ted. In brightness temperature, the two sigma cutoffs are 
0.20 K, 0.80 K, 0.14 K, and 0.26 K, for CO (3-2), CO (1- 
0), HCO+(4-3), and HCO+(1-0) respectively. Since the 
HCO + (4-3) emission was only detected in the central - 
100 to 200 km s _1 , only the densest gas in the central 
< 2" (1 kpc) is traced using this method. 

In order to understand the correlations seen in Fig- 
ure [9] better, a series of LVG calculations were per- 
formed. We adopt an LVG analysis that assumes a 
single cloud with spherical geometry, which requires 
three independent parameters as inputs: the H2 den- 
sity (nH 2 ), kinetic temperature (Tk), and the molecu- 
lar abundance per unit velocity grad ient (X/(dv/dr)). 
The c ollision rates were ob t ained from lFlower &; Launavl 
(fl985|) for CO and lFloweri (fl999f ) for HCO+. The quan- 
tity X/(dv/dr), which is different for different species, 
is physically the most ambiguous parameter. In or- 
der to decrease the number of unknown parameters, we 
adopt a constraint to the relative ab undance between 
CO a nd HCO+ (X HCO+ ~ 10- 4 X CO ; iBergin fc Laneren 
|1997|) obtained from chemical evolution models. By 
varying the CO abundance per unit velocity gradient 
from X CO /(dv/dr) = 10~ 4 to 10~ 9 (km s" 1 pc" 1 )- 1 , 
Xco/(dv/dr) ~ 10~ 7 (km s _1 pc -1 ) -1 provides the 
most sensible explanation of the data. Changing 
X C o/{dv/dr) to 10~ 6 (km s" 1 pc" 1 )" 1 will shift the 
temperature (vertical) contours leftward (see Figure [TO]) , 
yielding most of the molecular gas with temperatures 
in the range 100 - 200 K. This is higher than the cold 
and warm dust temp eratures pred i cted from a three 
component model by lArmus et al.l (|2006f ) where they 
find cold (27 K), warm (81 K) and hot (680 K) dust 
components to coexist over the entire NGC 6240 sys- 
tem. Using Xco I {dv / dr) = 10 -8 (km s _1 pc -1 ) -1 will 
shift the temperature contours rightward, yielding most 
of the data points in T < 40 K. This is lower than 
the values predicted in the nearby starburst NGC 253 
(| Jackson et aljfl995l ). but consistent with t he cold com- 
ponen t dust in NGC 6240 suggested by lArmus et al.1 
(2006), assuming that the gas and dust are in thermal 
equilibrium. 

Assuming that the gas temperature of 40 - 100 K is 
a reasonable representation of the molecular clouds in 
NGC 6240, then X C o/{dv/dr) ~ 10~ 7 (km s" 1 pc" 1 )" 1 
provides the most sensible abundance that explains 
the data in all channels with significant emission (Fig- 
ure [5]). A tight constraint to the H2 density of jjh 2 = 
2.Q5.0-5.4 cm~ 3 is obtained. The derived opacities range 
from t = 0.01 - 0.2 for the CO (1-0) line, r = 0.2 

2 for the CO (3-2) line, and t — 0.5 — 1 for the 
HCO + (4-3) line. The model predicts that the relative 
population of the HCO + (1-0) line be inverted, thus giv- 
ing negative opacities and excitation temperatures. The 
excitation temperatures (T ex ) derived from these analy- 



ses show that both CO lines are nearly thermalized (i.e. 

T ex ,CO(3-2) ~ T OXiCO (i-o))- 

Once these physical quantities are derived, it is now 
possible to constrain the beam filling factor of the CO 
emitting molecular gas. The observed line intensities are 
diluted by the ~ 2" (1 kpc) beam, and thus the true flux 
is obtained by scaling the observed flux by the reciprocal 
of the filling factor. Figure [Tl] and [12] show the correla- 
tion between T^o (3-2) and T^o (1-0) m each 21 km s _1 
channel where the emission of each line is significant (i.e. 
the central -133 to 187 km s _1 ). The filling factors of 
the CO lines, which are assumed to be the same between 
CO (3-2) and CO (1-0) as well as in all velocity channels, 
are varied from 1 (no scaling; Figure ITTj) to 0.2 (multiply 
the observed brightness temperatures by 5; Figure [T2")) . 
The lower limit to the filling factor of 0.2 is obtained 
by adjusting the filling factor until the highest observed 
points reach tt,h 2 ^ 10 54 cm -3 , a constraint given by the 
Rco - R-HCO+ pl°t m Figure [H A CO filling factor of 
unity gives the bulk of the gas in the temperature range 
of 40 - 200 K and density of n Ha ~ 10 4 cm" 3 . In this 
case, while the range of temperature is consistent with 
our earlier analysis (see above), the H2 density is an or- 
der of magnitude lower. A CO filling factor of 0.2 also 
yields a temperature range of 20 - 100 K, and the range 
of H2 density is also comparable to our earlier analysis. 
These results suggest that the true CO filling factor is 
between 0.2 to 1, but probably closer to 0.2. A low value 
of the filling fac tor would be inconsi s tent w ith the pre- 
vious studies by iDownes fc Solomon! (|1998[ ) where they 
suggest that the diffuse inter-clump medium may domi- 
nate the CO emission in ULIRGs, yielding a filling factor 
of unity. 

An estimate of the CO beam filling factor allows us 
to derive the volume filling factor of the C O emitting 
gas (see Equation 3 of lSakamoto et al]|1994f h Assuming 
that the CO (3-2) emission arises from GMCs of order 
50 pc each, and that the CO (3-2) emission region along 
the line of sight extends 4 kpc, the range of possible vol- 
ume filling factor is 0.004 - 0.05 for filling factors 0.2 - 1. 
Changing the size of the GMC to 25 pc (10 pc) gives 0.002 
- 0.03 (0.001 - 0.01) for the volume filling factor. Using 
the average density estimate from our LVG analysis (i.e. 
riH 2 ~ 10 5 cm -3 ), the derived molecular gas mass is 
then 10 9 " 10 M© in the central 0.5 kpc (r=0.25 kpc) re- 
gion. Our mass estimate is consistent with 6.9 x 10 9 M© 
derived using X=0.8 M (K km s^ 1 )" 1 pc 2 (see §3.1), 
and with (1 — 4) x 10 9 M© der ived in the central .5 kpc 
using CO (2-1) line emission (jTacconi et al.l ll999). As- 
suming that the average density in the larger spatial ex- 
tent is also 7Jh 2 = 10 5 cm" 3 gives ~ 10 10 M© in the 
central 1 kpc (r=0.5 kpc), and ~ 10 11 M© in the cen- 
tral 2 kpc (r=l kpc) region. The mass estimate given in 
a larger region (i.e. 2 kpc) is probably an overestimate 
since the molecular gas density in the periphery of the 
nucleus is likely lower than the central 0.5 kpc. Since 
the density derived from the LVG analysis only accounts 
for the densest gas in the central region, a proper mass 
estimate must consider the density variation of the cloud 
in the entire 1175 km s _1 velocity width, which is not 
possible here because of the varying mass sensitivities 
among the four data cubes. 

The LVG calculation also provides an estimate 
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of the H2 column density from the relation 3.0 x 
10 18 (X CO /(dv/dr)) n H2 = N co /dv. By adopt- 
ing X CO /{dv/dr) = 1(T 7 (km s' 1 pc' 1 )' 1 and 
n H2 = 10 5 cm" 3 , we find N co /dv = 3.0 x 
10 16 cm~ 2 [km s -1 ] -1 . Using the standard ab undance ob- 
serve d in galactic sources (Xco /Xh 2 = 10~ 4 ; Blak e et all 
I1987D and dv = FWZI = 1175 km s" 1 , we obtain Njj 2 = 
4.0 x 10 23 cm -2 . This is in excellent agreement with our 
previous column density estimate derived using the peak 
CO (3-2) flux. 

The peak brightness temperatures obtained in the 
CO (3-2) and 880^m observations (see §3.1 and §3.3) 
allow a confirmation of the properties derived above. In 
general, the brightness temperature and the true excita- 
tion temperature are related by Tb = f T ex (1 — e~ r ), 
where T B , /, T ex and r are the observed brightness 
temperature, beam filling factor, excitation temperature, 
and opacity respectively. Assuming the observed molec- 
ular clouds and the CO excitation are in thermal equilib- 
rium coupled with dust (i.e. T ex ~ T# ~ Tdust ~ 49 K), 
and using the range of filling factors of / = 0.2 - 1 (see 
above), we obtain r = 0.2 - 1.2, which is consistent with 
the opacities derived above. This is further evidence that 
the observed molecular clouds are moderately optically 
thin, and that they are close to thermal equilibrium. A 
similar analysis is performed on the dust continuum, and 
we derive the dust opacity Tdust — 0.003 - 0.01 using the 
same filling factors as CO, suggesting that the cold dust 
traced in 880^m is indeed optically thin. 

In summary, the important result obtained from these 
quantitative analyses is that the peak of the molecular 
gas in the central 1 kpc of NGC 6240 is warm, dense, and 
moderately optically thin. The derived temperatures and 
densities of NGC 6240 are consistent with the high values 
found in the nearby starburst galaxies M82 and NGC 253 
(IWild et al.lll992tlJackson et alJll995llSeaquist fc Fraverl 
2000). The low opacities predicted from these analyses 
are somewhat consistent with earlier results that inves- 
tigat e the opacities in M82 (r = 0.5 — 4.5; Ma o et ahl 
2000), but our new results predict much lower opaci- 
ties than previously sugg ested in ULIRGs (r = 3 — 10; 
iDownes fc Solomon! 11998 ). The low opacities found in 
the central region of NGC 6240 are possibly due to large 
velocity gradients caused by extremely turbulent gas, al- 
lowing emission in the core of the molecular cloud to 
escape efficiently. In addition, our predicted peak den- 
sities (nn 2 = 10 5 - 0-5 ' 4 cm" 3 ) are much higher than the 
range of values found in ULIRGs us ing low J transition 
CO e mission (nn 2 = 10 2 ' 3 ~ 4,3 cm~ 3 ; |Pownes fc Solomon] 
1998). The disagreement is likely due to the difference 
in the excitation conditions of the o b served molecular 
gas tracers; i.e. IDownes fc Solomonl (|1998f ) used cold 
(T ~ 10 K) and low density (n cr n ~ 10 3 cm' 3 ) trac- 
ers whereas our new observations trace much warmer 
(T - 30 K) and denser (n crit ~ 10 4 cm' 3 for CO (3- 
2)) molecular gas. This study, therefore, demonstrates 
the importance of using high density submm line trac- 
ers such as CO (3-2) and HCO + emission to study the 
physical conditions in the nuclear ISM in ULIRGs. 

4.2. The Gas to Dust Mass Ratio 

The gas-to-dust mass ratio (i.e. M^/Mj) provides 
an important measure of the relative abundance be- 
tween gas and metallicity. An average Mn 2 /Md over 



the entire galaxy is often derived in single dish work, 
where i t is found that M h^/M^ ~ 100 for Galactic 
sources (|Hildebrandlll983l). MuJMd ~ 200 - 300 for lo- 
cal L IRGs /ULIRGs jContini fc Continil l200l lYao et all 
120031 ISeaauist et ap|2004j), and M H ,/Md = 15 - 231 in 
high-z sources (jSolomon fc vanden Bout! [20051 ). Wilson 
et al. (2006 in prep) finds M H2 /M d = 357 ± 95 from 
a sample of 13 LIRGs/ULIRGs observed in high resolu- 
tion with the SMA. It is found that the global gas-to-dust 
mass ratio in NGC 6240 is M H ,/M d = (147 ± 57) (see 
§3.3). 

Our new high angular resolution CO (3-2) and 880 fim 
maps allow us to study the spatial variation of MhJM^. 
We have converted the CO (3-2) and 880/im maps to 
Mh 2 and Md maps. The gas-to-dust mass ratio map is 
shown in Figure [13] The peak of Mh 2 /M<z occurs be- 
tween the two AGNs, and near S2, suggesting an over- 
abundance of dense gas relative to dust in these regions, 
assuming that most of the 880^im emission is tracing 
dust. The high ratios seen north of the northern AGN 
are due to low dust masses at the edge of the 3<r cutoff. 
The peak ratio is ~ 150. We note that the ratio map ob- 
tained here is only available in the central kpc where the 
dust emission is significant (> 3a). The CO (3-2) emis- 
sion is much more extended than the dust emission and 
hence the average Mh 2 /M<2 in Figure [13] is lower than 
the global Mh 2 /M,i (147 ± 57). Higher ratios are seen 
in the north-south direction (i.e. Mn 2 /Md ~ 100 - 160) 
along the two AGNs than the east-west direction (i.e. 
Mh 2 /M,j ~ 30 - 100). The steep decline in the wings of 
dust emission may suggest substantial free-free emission 
associated with the starburst winds, or a local gradient 
in the metallicity. 

4.3. The FIR to CO (3-2) Luminosity Ratio 

It is empirically known that the ratio between FIR and 
CO (1-0) luminosit ies is higher for ULIRGs than for nor- 
mal sp iral galaxies (|Solomon et al.|[l997[) . ISolomon et al.l 
(|1997| ) attribute this result to the presence of com- 
pact, optically thick dust emission, which in turn yields 
higher dust temperatures than in the optically thin limit. 
However, these measurements are performed using the 
CO (1-0) transition, and it is found that this value 
is different when the ratio between FIR and CO (3- 
2) is taken (i.e. compare log Lp^/log Lc Q (3-2) — 
1.43 ± 0.08 wit h log L Fm /log L co( i_ ) = 1-72 ± 0.21; 
lYao et al.ll2003l) . The FIR-to-CO (3-2) luminosity ra- 
tio appears to be even lower in high-z sources where 
the FIR luminosity is typically an order of magni- 
tude larger (i.e. Lfir ~ 10 13 L Q ). From the high- 
er sources with a robust CO (3-2) detection, the ratio 
is found to be log Lpm/log L cq( 3 _ 2 ) = 1.23 ± 0.03 
(|Solomon fc vanden Boutl l2005t llono et alj I2006D . Us- 
ing the CO (3-2) integrated intensity from our full res- 
olution map, and the IRAS measurement for the FIR 
luminosity, the ratio for NGC 6240 is found to be 
log L Fm /log L CO (3-2) = 1-16. 

NGC 6240 thus contains much more molecular gas 
relative to the intensity of starburst activity than sim- 
ilarly FIR bright systems observed in the local universe. 
The ratio is even lower than the values derived in high- 
z sources. The overwhelming amount of molecular gas 
available in the central region of NGC 6240 suggests that 
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there is enough gas to support star formation for at least 
another ~ 10 8 years, assuming that all of the gas will 
eventually be converted to stars with a constant star 
formation rate of 60 M yr _1 . The wealth of molec- 
ular gas may also suggest that the star formation in 
NGC 6240 will increase to a level comparable to that 
seen in Arp 220. 

5. DISCUSSION 

5.1. The Nature of the Western Molecular Complex 

(WC) 

Recent galaxy-ga laxy collision simulations (e.g. 
Hopk ins et al.l l2005f ) that trace the evolution of the 
central AGN have shown that significant outflow ac- 
tivity may be obse rvable in high J CO transitions 
up to J = 6 - 5 (|Naravanan et al.l I2006D . While a 
massive (3.5 x 10 s M ) off-nuclear molecular complex 
such as WC is often explained by starburst superwinds 
(jHeckman. Armus. fc Milevl Il990h . we investigate here 
the possibility that the origin of WC is gas entrained in 
outflow from the AGN of one of the progenitor galaxies. 
From an estimate of the cloud diameter of D = 2.7 kpc 
from the CO (3-2) map and the velocity dispersion of 
a v = 400 ± 100 km s _1 estimated from the FWHM 
of the CO (3-2) line in WC, we find a virial mass of 
M V i r = 3.0 x 10 11 M . The error in this virial mass 
estimate could be as large as a factor of (2-3) due 
to large uncertainties in the cloud size and velocity 
dispersion. By assuming the Galactic dynamical mas s 
(M dyn ) to M Ha ratio of 136 (|Nakanishi fc Sofud l2006f ) . 
we estimate a total mass of M tot = 4.8 x 10 10 M for 
WC. Although the uncertainties in these measurements 
are large, the derived virial mass is an order of magni- 
tude larger than the total mass of WC. This suggests 
that WC is not bound by its own gravitational force and 
will likely disperse within D/a v ~ 7 x 10 6 years. The 
systemic velocity of WC is v sys = 7449 km s" 1 , which 
yields a relative velocity wvvc ~ 210 km s~ x assuming 
that WC is ejected from the northern nucleus. Using 
a projected distance between the northern nucleus and 
WC of 4.4 kpc, the timescale for WC to traverse to the 
current location is ~ 2 x 10 7 years. 

Although all of these timescales contain significant un- 
certainties, the short dispersion time (7 x 10 6 years) sug- 
gests that it is unlikely that the cloud would have sur- 
vived the estimated transversal time of 2 x 10 7 years 
to reach the current locat ion. Simulations conducted 
by iNaravanan et al.1 (|2006f ) predict that the total gas 
mass of an ejected cloud could be as large as ~ 10 s M 
with a cloud disp ersion time of t ~ 2.9 x 10 8 years 
(|Naravanan et al.l 12006). While the CO (3-2) derived 
gas mass is consistent with the ejected gas mass pre- 
dicted from simulations, the derived dispersion timescale 
is a factor of 40 shorter. These results may suggest that 
the AGN outflow scenario is less appealing to explain 
WC, and other scenarios such as molecular outflow from 
starburst superwinds may be more appealing to explain 
the large and extended molecular gas. The starburst 
outflow scenario to explain the western extension of ra- 
dio continuum emis sion in NGC 6240 was suggested by 
iColbert et all (fl99l . 

Possible scenarios to explain WC are not limited to 
starburst superwinds or to AGN outflows and both of the 



following are also possibilities; (1) molecular gas ejected 
from the violent arm-arm collision, (2) molecular gas 
complex formed locally as a result of gravitational insta- 
bility. Ejection of molecular gas is possible as arm-arm 
collision could occur in each pericentric passage. How- 
ever, ejecting a molecular gas mass of order 10 8 M and 
diameter 1-2 kpc from the potential of a tidal arm may 
require a significant amount of energy. While it is pos- 
sible that gravitational collapse could form a CO cloud 
from an ambient H I cloud, it is already demonstrated 
that the virial mass is an order of magnitude larger than 
the estimated total mass of WC. Hence these two scenar- 
ios are less appealing to explain the origin of WC than the 
starburst superwind or the AGN outflow scenarios. Our 
current data, however, is not sensitive enough to address 
the exact origin of WC, and future observations with 
higher angular resolution and sensitivity should provide 
us further understandings to the possibilities of these sce- 
narios. 

5.2. The Central 1 kpc of NGC 6240 

Our new analysis presented in §4.1 shows that the 
molecular gas in the central 1 kpc of NGC 6240 is warm, 
dense and moderately optically thin. These results sug- 
gest that the excitation conditions of molecular gas sur- 
rounding the two AGNs is similar to the nuclear regions 
of nearby starbursting galaxies. The low opacity is likely 
caused by highly turbulent gas, which is also supported 
by the dominance of non-circular kinematics entangled 
with possible disk-like rotation around the northern nu- 
cleus (see §3.1.2). The southern nucleus shows little evi- 
dence of disk-like rotation possibly due to the dominance 
of non-circular motion, or due to a face-on orientation of 
a rotating disk. Although the projected nuclear sepa- 
ration is different, the double nucleus configuration is 
similar to that seen in Arp 220 where a merger of two 
counter-rotating disks is suggested from high resolution 
CO (2-1) observations ([Sakamoto et al.l 1999). Higher 
angular resolution observations should provide further 
insight to these possibilities. 

In addition to the large amount of dense molecular 
gas, there is significant multi-wavelength evidence that 
the central 1 kpc also harbors substantial amounts of 
dust. While we caut ion here that emiss ion traced in 
880/im and 1.3 mm (jTacconi et al.l [T999[ ) are possibly 
more dominated by non-thermal emission, the ubiqui- 
tous presence of the optically thin dust emission in the 
medium between the two AGNs is a significant evidence 
of cold dust. Further, from comparison between star for- 
mation activity traced in optic a l and the AGN activity 
traced in X-ray, iGerssen et al.l (|2004l ) have shown that 
the two nuclear star forming regions in both progenitor 
disks are not coincident with the AGN positions but are 
offset by ~ C/5, suggesting that foreground dust extinc- 
tion in the central region is significant. The dusty central 
1 kpc also shows significant evidence of a coexistence of 
hot and shocked gas as seen in detections o f maser emis- 
sion fHagiwara. Diamond & Miyoshi 2003) and H2 emis- 
sion (jOhvama et al.l [2000h . Both of these are detected 
near SI and S2, providing evidence of shocked gas asso- 
ciated with either the AGN jet or the nuclear starburst 
activity. Summarizing the multi- wavelength properties, 
the central 1 kpc of NGC 6240 harbors two AGNs sur- 
rounded by a wealth of warm and dense molecular gas, 



as well as hot (T = lO^ 3-7 -* K) gas and a large amount 
of dust, and the region is undergoing a high rate of star 
formation activity. 

6. SUMMARY 

We present high resolution CO (3-2), HCO+(4-3) 
and 880/xm imaging of the luminous infrared galaxy 
NGC 6240 obtained at the SMA. We provide our main 
findings below; 

1. The distribution of CO (3-2) is extended on a 4 kpc 
scale, but the HCO + (4-3) emission is more com- 
pact and concentrated in the central 1 kpc between 
the two nuclear AGNs. The CO (3-2), HCO+(4- 
3), and 880^m emission all peak between the two 
AGNs where the H 2 emission was detected previ- 
ously. 

2. The kinematical information provided by the 
CO (3-2) emission shows a rotating disk-like fea- 
ture centered around the northern AGN, but the 
kinematics of gas between the two nuclei is ex- 
tremely turbulent. 

3. Our LVG analysis shows that the gas near the 



peak of the continuum emission is warm (T = 20 
- 100 K), dense (n H2 = lO 50 " 5 - 4 cm" 3 ) and mod- 
erately optically thin (r = 0.2 - 2) in the central 
1 kpc of NGC 6240. 

4. We propose that the isolated molecular emission 
with 3.5 x 10 8 M Q seen west of the main CO 
component may be associated with outflows from 
starburst superwinds, but the gas outflow scenario 
from one of the central AGN is not completely ruled 
out. 
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TABLE 1 
Derived Properties 



Observation 


Property 


Value 


CO (3— 2) 








Vsys (km s L ) 


7449 ± 25 




FWZI (km s _1 ) 


1175 ± 25 




Sj,oV (Jy km s _1 ) 


2834 ± 20 




peak R.A. 


16 h 52 m 58.9 s 




peak Deck 


2°24'03'.'7 




peak T B (K) 


8.6 


HCO+(4-3) 








Vsys (km s L ) 


7339 ± 50 




FWZI 1 (km s- 1 ) 


300 ± 50 




Svdv 1 (Jy km s- 1 ) 


28 ±4 




peak R.A. 


16 h 52 m 58.9 s 




peak Deck 


2°24'03'.'8 




peak T B (K) 


0.5 


880/im continuum 








S 8 so (mJy) 


105 ±6 




peak R.A. 


16 h 52 m 58.9 s 




peak Deck 


2°24'03'.'4 




peak T B (K) 


0.1 



1 The FWZI and Svdv of the HCO+(4-3) emission may be 
underestimated since the emission occurs near the edge of 
the bandpass (sec text). 
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Right Ascension (J2000) 

Fig. 1 — (left) Digital Sky Survey image of NGC 6240 overlaid with the 35" SMA field of view, (right) The HST NICMOS 2.2/xm image 
of the central 2" (1 kpc) of NGC 6240 in grayscale overlaid with various peak positions from radio/mm/submm observations. Since accurate 
astrometry is not provided in the HST image, we align the peak of the southern 2.2/^m nucleus with the 2.4 GHz nucleus. The length of 
the crosses indicate the astrometric u ncertainties of each m e asure ment. The points labele d a and b are t he loc ations of the two AGNs. 
Coord inates are obtained from (a,b,c) Gallimorc & Bcswick (2004), (d,f,h) this work, (e,g) Tacconi ct al. (1999), and (i) Nakanishi et al. 
120051) . 




15 10 5 -5 -10 -15 1 15 10 5 -5 -10 -15 

RA offset (arcsec; J3000) RA offset (arcsec; J2000) 



Fig. 2. — The synthesized beam using (left) natural and (right) uniform weighting. The contour levels are 10 - 90% of the peak. 
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10 5 -5 -10 10 5 -5 -10 10 5 -5 -10 10 5 -5 -10 
RA offset (arcseo; J2000) RA offset (arcsec: J2000) RA offset (arcsec: J8000) RA offset (arcsec: J2000) 

Fig. 3. — The natural weighted CO (3-2) channel maps in the velocity range 6889 to 8014 km s . Although the original channel 
maps were made with 25 km s — 1 velocity resolution, the maps presented here are made by averaging three channels (75 km s _1 ) in order 
to present the extremely broad velocity width in a single figure. The velocities are shown in the upper left corner of each panel, and the 
synthesized beam is shown in the lower left. The two crosses show the location of the two AGNs (Nl and SI in Figure [TJ. The contour 
levels are 16 mjy (1 sigma) (x -5,-3,3,5,7,9,15,20,40,60,80,100,120,140). 





16 h 52 m 59?10 59?05 59?00 58?95 58?90 58?85 58?80 58?75 
RA (J2000) 

Fig. 4.— (left) The CO (3-2) moment map of the central 18" of NGC 6240. The contour levels are 5.5 Jy km s" 1 (x 
3,5,7,9,15,20,30,50,100,200). We define the isolated CO (3—2) emission complex west of the main CO (3—2) concentration as WC. (right) 
The uniformly weighted CO (3-2) moment map of the central 6" of NGC 6240. The contour levels are 16 Jy km s' 1 (x 3,5,7,9,15,21,27). 
The crosses are the locations of Nl and SI (see Figure[TJ, and the synthesized beams are shown in the lower left. 
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i 52 m 59 3 10 59 a 05 59 a 00 5B S 95 58 a 90 58 a B5 5B a B0 5B?75 
RA (J2000) 



Fig. 5. — {left) The CO (3-2) velocity distribution (moment 1) map from the natural weighted map in the central 14" of NGC 6240. 
The velocity contours are spaced at 50 km s — 1 . (right) same as left but using uniform weighting. The corresponding moment maps are 
shown in Figure [4] The units are in km s _1 . 
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16 h 52 m 58?6 



58 B 5 



58*4 



RA (J2000) 

Fig. 6. — The CO (3—2) velocity distribution of the western complex (WC) alone. The velocity contours are spaced at 30 km s~ 1 . The 
units are in km s _1 . 
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-4-2034 -2-1012 -2-101 

POSITION POSITION POSITION 



Fig. 7. — The position velocity diagram for (left) the 2'.'0 X l'.'S resolution map (natural weighting). The x-axis shows the offset position 
from the center in arcseconds, and the y-axis shows the velocity in km s -1 . The contours represent 10,20,30,40,50,60,70,80,90% of the peak. 
These PVDs are made by slicing the cube along the two AGNs (Nl and SI) and centered exactly halfway between the two AGNs. The 
location of the two AGNs are shown as vertical lines. The linewidth near the two nuclei is extremely large, with emission slightly lopsided 
toward the redshifted velocity, (middle) Similar to left but for the l'.'l X 0'.'8 resolution map (uniform weighting). The contours represent 
the 1,5,10,20,40,60,80% of the peak, and (right) 0'.'9 X 0''5 resolution map (using the longest baselines only). The contours represent the 
15,30,45,60,75,90% of the peak. The nuclear region is resolved into two distinct components in the highest resolution map. 




16 h 52 m 59?l 59'0 5S 3 9 58'8 56?7 16 h 52 m 59 a l 59 s 58*9 5S'B 5B'7 



RA (J2000) RA (J3000) 

Fig. 8. — (left) The integrated intensity map of the HCO"'"(4-3) emission. The contours are 5 Jy km s — 1 (x 3,4,5,6,7,8). (right) The 
880^tm emission map. The contours are 6 mjy (x -4,-3,3,4,5,6,7). Both of these maps are made using natural weighting, and the symbols 
in both maps are the same as in Figure [4] 
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Fig. 9.— The R C o(CO (3-2)/CO (1-0)) vs. R HC o+ (HCO+(4-3)/HCO+(l-0)) plot for each 56 km s^ 1 channel. The crosses show 
the pixel-to-pixel correlation of the two ratios. The cross shown in the lower right of the first panel represents a typical error bar associated 
with the data. Xco/(dv/dr) = 10~ 7 (km s _1 pc -1 ) -1 is adopted for this particular model. The solid line shows the temperature contours 
of 40, 60, 100, and 200K, and the dashed line shows density contours of 10 5 and 10 5 4 cm -3 . The lower density contours (i.e. 10 4 cm -3 ) 
lie in the lower left corner of each panel, which are not visible with the current scaling. 
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Xco/dv/dr = 10 Xco/dv/dr = 1 Xco/dv/dr = 1 




R CO 



Fig. 10.— The model R C o(CO (3-2)/CO (1-0)) vs. R HC o+ (HCO+(4-3)/HCO+(l-0)) shown for different X co /(dv/dr) for compari- 
son. The crosses show the data points from v = 7298 km s _1 in Figure[5] The solid line shows the temperature contours of 40, 60, 100, and 
200K, and the dashed line shows density contours of 10 5 and 10° A cm -3 (and 10 4 5 cm -3 for X t 7o/(dv/dr) = 10 -6 (km s _1 pc -1 ) -1 ). 
Using Xco/(dv/dr) = 10 — 6 (km s — 1 pc -1 ) -1 will shift the temperature (vertical) contours leftward, yielding most of the molecular gas 
with temperatures in excess of 200 K. Using Xpo/^dv/dr) = 10 — 8 (km s _1 pc -1 )" 1 will shift the temperature contours rightward, yielding 
a wide range of temperatures but mostly T < 60 K. 




Fig. 11. — The pixel-to-pixel correlation between Tco(3-2) an d Tco(l-O) ( m Kelvin) are shown in crosses, where a beam filling factor 
of unity is assumed. The solid lines show the T^ = 20, 40, 60, 100, 200 K, and the dashed line shows nn 2 = 10 4 cm -3 contour. 
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Fig. 12 — 



1 CO(l-O) 

Similar to Figure [TT1 but a beam filling factor of 0.2 is assumed. 
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